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Microscopic dynamics of the oblate-prolate shape coexistence/mixing phenomena in ^*Se 
and ^^Kr are studied by means of the adiabatic self-consistent collective coordinate (ASCC) 
method in conjunction with the pairing-plus-quadrupole (P+Q) Hamiltonian including the 
' quadrupole pairing interaction. Quantum collective Hamiltonian is constructed, and exci- 

^ ' tation spectra, spectroscopic quadrupole moments and quadrupole transition properties are 

^ , evaluated. The effect of the time-odd pair field on the collective mass (inertia function) of 

the large- amplitude vibration and the rotational moments of inertia about three principal 
axes is evaluated. Basic properties of the shape coexistence/mixing are well reproduced. 
, The calculation indicates that the oblate-prolate shape mixing decreases as the angular 

Jf\ . momentum increases. 

a^ 

^ ' §1. Introduction 

Microscopic understanding of nuclear collective dynamics is one of the goals of 
' nuclear structure theory. The quasiparticle random phase approximation (QRPA) 

■ based on the Hartree-Fock-Bogoliubov (HFB) mean-field is the well known the- 
oretical approach to the collective dynamics, but it is applicable only to small- 
amplitude collective motions around the local minima of the potential energy sui- 
face.'^''^''^®''^ Nuclei exhibit a variety of large-amplitude collective process such 

5^ \ as anharmonic vibrations, shape coexistence, and fission. Though the challenge 

to construct microscopic theory of large-amplitude collective motion has long his- 

t ory,El' 'O M Mi ,llD ,111 ,113 ,111 ,113 ^ Mi Mi Mi Mi ,[1D ,123 ,|23]i ^ Mi ,123 ,|27li ,|28]i Mi Mi .ED .133 ,|33]i ,[31 ,|3g S Mi S ,E 
some serious problems remain unsolved. 

The self-consistent collective coordinate (SCO) methocP^''^ is a microscopic 
theory of large-amplitude collective motion. This method with the {"/]* ,rj) expansion 
technique enables us to extract the collective variables from the huge-dimensional 
phase space associated with the time-dependent Hartree-Fock-Bogoliubov (TD- 
HFB) state vectors, and to derive the collective Hamiltonian starting from a micro- 
scopic Hamiltonian. The SCO method has been successfully applied to various kinds 
of non-linear phenomena in nuclei; such as anharmonic 7-vibrations,H3''^'H3'HHli 
shape phase transitions,S3''^'H21'''53 two-phonon state^SD and collective rotations. 'ISSS 
The (??*,??) expansion is not suitable, however, for large-amplitude collective mo- 
tions like shape coexistence/mixing phenomena, where microscopic description of 
many-body tunneling effect between different local minima in the collective poten- 
tial energy surface becomes the major task. 
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The adiabatic SCC (ASCC) method^^^ is an alternative way to solve the basic 
equations of the SCC method assuming that the large-amplitude collective motion of 
interest is slow (adiabatic). Under this assumption, the basic equations of the SCC 
method are expanded up to second order with respect to the collective momentum, 
but the collective coordinate is treated non-perturbatively. Quite recently, we have 
given a rigorous formulation of the ASCC method in which the gauge invariance 
with respect to the particle number fluctuation degrees of freedom is taken into 
account. 

As is well known, in contrast to the remarkable progress in the calculation of 
the collective potential energy, the present status of the microscopic theory is quite 
unsatisfactory concerning the evaluation of the collective mass (inertia function) 
associated with the collective kinetic energy. Although the Inglis-Belyaev cranking 
mass is widely used, it violates the self-consistency by ignoring the effect of the 
time-odd component of the moving mean-field.'^ The time-odd mean-field effect is 
taken into account in the collective mass derived by the QRPA, but its application 
is restricted to small-amplitude collective motions around the equilibrium states. 
Concerning large- amplitude collective motions, though the effect of the time-odd 
component generated by the residual particle-hole interaction was investigated in a 
few decades ago,'^ the time-odd effect generated by the residual pairing interaction 
has not been discussed so far. Quite recently, we have shown, using the ASCC 
method in conjunction with the schematic model Hamiltonian,'^''^'!^''® that the 
time-odd pair field increases the collective massW^ It remains to be seen, however, 
how it affects the shape coexistence dynamics discussed below. 

Let us turn to the recent experimental data we are interested in. Nuclei along 
the N = Z line evolves its shape drastically with changing number of protons and 
neutrons.lSS'El'llS'lSB'ES The HFB calculatioiJSB indicates that various shapes wih 
appear along the N = Z lines: a triaxial ground state for ®^Ge, oblate ground 
states for ^^Se and ^-^Kr, strongly deformed prolate ground states for ''^Sr, ^''Zr and 
^^Mo. Furthermore, oblate and prolate states may coexist in these nuclei except 
^^Ge. In ^^Se and ^^Kr, the ground and excited states corresponding to the oblate 
and prolate shapes have been found experimentally.'S3'l53'[S3 Prom the viewpoint of 
collective dynamics based on the mean-field theory, it is expected that the oblate 
and prolate shapes are mixed by the many-body tunneling effect through the poten- 
tial barrier lying between the two local minima in the potential energy landscape. 
The low-lying states in ^'^Se and ^■^Kr have been investigated by various theoretical 
approaches beyond the mean-field approximation: Large-scale shell model calcula- 
tion for ®^Se using the pfg-shells outside the ^^Ni cove^^ shell model Monte Carlo 
calculation for ''^Kr employing the pf-sdg shells,"^ configuration mixing calculation 
for 72-78j^j, ^j^g basis of the particle number and angular momentum projected 
generator coordinate method,^-^ and Excited Vampir variational calculation for ®^Se 
and 72Kr.llD.El.E3 Quite recently, Almehed and WalelP'E3.l33.[3i discussed collec- 
tive paths connecting the oblate and prolate minima in ^^Se and 72-78j^j. j-^y jng^j^s 
of the approach similar to the ASCC method. 

The ASCC method was first testecP^ in the schematic modePSJ '121.123 .ED and 
then appliecpl to the oblate-prolate shape coexistence phenomena in ^^Se and ^^Kr 
with the use of the pairing-plus-quadrupole (P+Q) Hamiltonian.E3.E3.E3.E3 Jn 
both nuclei, the one-dimensional collective path connecting the two potential local 
minima is extracted. It was found that the collective path runs approximately 
along the valley of the potential energy surface lying in the triaxial deformed region. 
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This indicates that the triaxial degree of freedom is essential for the description of 
large-amphtude shape mixing in ^'^Se and ^^Kr. In the previous workf^ however, 
requantization of the collective Hamiltonian was not done, and excitation spectra, 
electromagnetic transition probabilities, and shape mixing probabilities in individual 
eigen-states were not evaluated. 

This paper presents the result of the first application of the gauge invariant for- 
mulation!^ of the ASCC method to nuclear structure phenomena. Thus, the major 
effort is directed to examine the feasibility of the gauge-invariant ASCC method for 
describing the shape coexistence/mixing phenomena. Hereafter we call this new ver- 
sion "ASCC method" dropping the adjective " gauge-invariant" for simplicity. More 
detailed investigation of experimental data and comparison with other approaches 
are planned for future. We derive the quantum collective Hamiltonian that de- 
scribes the coupled collective motion of the large-amplitude vibration responsible 
for the oblate-prolate shape mixing and the three-dimensional rotation of the triax- 
ial shape. To evaluate the rotational moments of inertia, we extend the well-known 
QRPA equation for rotational motion, which yields the Thouless-Valatin moment 
of inertia, to the non-equibrium states that are defined in the moving-frame as- 
sociated with the large-amplitude vibrational motion. To clarify the role of the 
time-odd pair field in shape mixing dynamics, we investigate, with the use of the 
P+Q Hamiltonian including the quadrupole pairing interaction, its effects on the 
collective mass of large-amplitude vibration, the rotational moments of inertia, the 
energy spectra, transition probabilities, and shape mixing probabilities in individual 
eigen-states. 

This paper is organized as follows: The basic equations of the ASCC method are 
summarized in ^ Quasiparticle representation of the microscopic Hamiltonian is 
given in ^ Solving procedure of the ASCC equations is presented in ^ Collective 
Schrodinger equation is derived in ^ Results of numerical calculation for energy 
spectra, spectroscopic quadrupole moments and quadrupole transition probabilities 
of low-lying states in ^®Se and ^^Kr are presented and discussed in ^ Concluding 
remarks are given in ^ 

A preliminary version of this work was previously reported in Ref. I8ip . 

§2. The ASCC method 

2.1. Basic equations of the ASCC method 

We first recapitulate the basic equations of the ASCC method. The TDHFB 
state \(l){t)) is written in terms of the collective variables as 

m)) = m,P,H^,n)) = e-E.^M^M m,P,n)) , (2-1) 

where q and p represent the one-dimensional collective coordinate and collective 
momentum, respectively. The variables ip = {(p^"-\ cp^P^) and n = {n^'^\n^^) denote 
the gauge angles in particle number space and number fluctuations, respectively, 
which correspond to the canonical coordinates and momenta of the pairing rotation 
restoring the particle number conservation broken by the HFB approximation. The 
operator iV^"^^ = iV^"^^ — n'^'' represents the particle number measured from the 

reference value N^^^ , which is set to the number of the valence protons (r = p) and 
neutrons (r = n) in the model space. 

The intrinsic state with respect to the pairing rotation is written as n)) = 
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^zG(q,p,n) where = \4>{QtP = 0, n = 0)). Assuming that large-amphtude 

coUective motion is adiabatic, that is, the collective momentum p and the number 
fluctuation n are small, we expand the one-body operator G{q,p,n) with respect to 
p and n^"^) and consider only the first order: 

mq,P,n)) = e^pQW+«E.nMe(^)(,) |^(^)^ ^ (2-2) 

where Q{q) is a time-even one-body operator, while 0^'^\q) is a time-odd one-body 
operator. Using quasiparticle creation and annihilation operators, aa{q) and aa{q), 
defined with respect to a moving-frame HFB state \<J){q)), which satisfy the condition 
d-aiQ) \4'{q)) = 0) they are written as 

Q{q) = Q%) + Q''{q) 

= E (Qa/3(9)4(9)4(^) + Qa*^il)MQMQ) + QU^)aUq)MQ)) , (2-3) 

a/3 

0(-)[q) = ^ [0^^l^{q)aUq)al{q) + 0S''*(g)a^(g)a«(g)) . (2-4) 

0/3 

Note that the operator Q{q) contains the i3-part (third term in the r.h.s.), in addition 
to the A-part (first and second terms) in order to satisfy the gauge-invariance of the 
ASCC equations.'^ In the following we will omit the index q in the quasiparticle 
operators for simplicity of notation. 

The basic equations of the SCC method consists of the canonical variable con- 
ditions, the moving-frame HFB equation, and the moving-frame QRPA equations. 
Below we summarize the lowest order expressions of these equations with respect to 
the expansion in p and n (See Ref. |56|) for their derivations). The canonical variable 
conditions are given by 



{<P{q)\P{q)\Hq))-- 


= 0, 


(2-5) 


{Hq)\Q{q)\m)-- 


= 0, 


(2-6) 


{cp{q)\N^^UHq))-- 


= 0, 


(2-7) 


(</>(g)|0M(g) 10(g)) = 


= 0, 


(2-8) 


(0(g)|[0W(g),7V(-')] 10(g)) = 




(2-9) 


(0(g)|[Q(g),0W(g)] 10(g)) = 


= 0, 


(2-10) 


imi^lHq))- 


= -1, 


(2.11) 



where P{q) is the local shift operator defined by 

P{q)\m)=^^^\m)- (2-12) 

Differentiating (j2-6p and (j2-7p with respect to q and using (|2-lip and (j2T2p . we 
obtain 



{cl){q)\[Q{q),P{q)mq))=i, 
(0(g)|[ivW,P(g)] 10(g)) =0. 



(2-13) 
(2-14) 
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Equations (12-51) , (I2-6P , (j2-7p and (|2-8p assure that the constant term of those oper- 
ators are zero in their quasiparticle representations, while Eqs. (j2-9p . (j2-10p . (j2-13p 
and (|2-14p guarantee orthonormahzation of the collective mode and the number 
fluctuation modes. Equation (j2-lip defines the scaling of the collective coordinate. 
The moving-frame HFB equation is given by 

5{m\HM{q)\m)=^, (2-15) 

where 

HM{q) = H-Y, A(^)(g)iVW - -^Q{q) (2-16) 

T 

represents the moving-frame Hamiltonian with the chemical potential )S^\q) and 
the collective potential V[q) defined by 

\^^\q)=^. ^ . .= {m\[H,^e^^\q)mq)). (2-17) 
dn^^i p=o,n=o,i=o 

Viq)=n{q,p,n,I) = {cl>{q)\ H {cj^iq)) . (2-18) 

p=0, 71=0,7=0 

The moving-frame QRPA equations are given by 

6 {4>{q)\ [HM{q),iQ{q)] - B{q)P{q) \<i,{q)) = 0, (2-19) 



5 {m\[HM{q),P{q)] - iC{q)Q{q) 

- ^^\[HM{ql^-^Q{q)],iQ{^^^ (2-20) 

where B{q) and C{q) represent the inverse collective mass and the local stiffness, 
respectively. They are defined by 

^(9) =0 !,=o,r^=0,r=0 = ^"^^^^l [[H,iQ{q)liQ{q)] \m) , (2-21) 

^/ , 5V I dBdV , , 

Note that the ASCC equations, (|2T5p . ()2-19p and (j2-20p . are invariant against 
the following transformation!^ 

Q{q)^Q{q)+a^^^N^^\ 
fIV 

AM(5)^AM(<i)-aW^(<,), 

(2.23, 

Therefore it is necessary to fix the particle number gauge for neutrons and protons 
in order to find the unique solution of the ASCC equations. The algorithm to 
find simultaneous solutions of Eqs. (j2-15p . (j2-19p and (j2-20p fulfilling the canonical 
variables conditions and the gauge-fixing condition is described in 21 
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In this paper, we take into account the rotational motion as weh as the large- 
amphtude vibrational motion by considering the collective Hamiltonian defined as 
follows: 

n{q,p,n,I) ={(j){q,p,n)\H\cl}{q,p,n)) + Yl^jr^Ii 



-V{q) + lB{q)p' + A(^)(g)nM + J] j^jf- (2-24) 



The first and the second terms represent the potential and kinetic energies of the 
large-amplitude collective vibration, respectively, while the third and the fourth 
terms the energies associated with the particle-number fluctuations and the three- 
dimensional rotations of triaxially deformed mean fields, respectively. The three 
rotational moments of inertia J^i{q) are defined with respect to the principal axes 
associated with the moving-frame HFB state \(f>{q)) and evaluated by 

6{(t>{q)\ [HM{q)A{q)] - \j^-\q)h \m) = 0, (2-25) 
{mm{Q)J3\\<^{<l))=^i3^ (2-26) 

where ^i{q) and /j represent the rotational angle and the angular momentum op- 
erators, respectively. These equations reduce to the well-known QRPA equations 
giving the Thouless-Valatin moments of inertid^ when \(t){q)) is an equilibrium 
state corresponding to a local minimum of the collective potential energy V{q). We 
call them "Thouless-Valatin equations," although they are in fact extensions of the 
QRPA equations for collective rotation to non-equilibrium HFB states \(j){q)). Note 
that Hm{q) appears in Eq. (j2-25p instead of H. 

§3. Hamiltonian 

We adopt the Hamiltonian consisting of the spherical single-particle energy, the 
monopole and the quadrupole pairing interactions, and the quadrupole particle-hole 
interaction: 



T K=-2 K=-2 



where the monopole pairing operator A^'^'^^ the quadrupole pairing operator B^^\ 
the quadrupole particle-hole operator D2K are defined by 

i«t = 44, (3-2) 

(fc,fe)6T 

^J? = E ^2^(^044' (3-3) 

Her 

^2/^= E E^22(^o4c^- (3-4) 

T=n,p kler 
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Here c\ is the nucleon creation operator, and k denotes the set of the quantum 

numbers of single-particle state {Nj., jk,lf^,mk)- The operator represents its time- 
reversal state, 

4 = (-iy'=+'"'=ct_„ (3-5) 

where the index —k denotes {Nj., j^^lj^^ —mi.). The quadrupole matrix elements are 
given by 

D^^j>{kl)=al{k\r^Y2K\l), (klGr), (3-6) 

where the factors, = {2N/A)'^/^ and = (2Z/y4)^/^, are multiphed to yield the 
same root mean square radius for neutrons and protons. For N=Z nuclei like ^^Se 
and ^■^Kr, these factors are unity. Following Baranger and Kumar,!^ we employ a 
model space consisting of two major oscillator shells (the total quantum number of 
the lower shell is denoted Nl and that of the upper shell Nl + 1), and multiply 
the reduction factor C = [N^ + 3/2) /{N^ + 5/2) to the quadrupole matrix elements 
D^]l{kl) of the upper shell. Following the conventional prescription of the P-l-Q 
model, we ignore the Fock terms. Accordingly, we use an abbreviation HB in place 
of HFB in the following. 

We rewrite Hamiltonian (j3-ip into the following form 



^ = E ^^4-^ - I E -sF^'-^F^^^ + \ , (3-7) 

k s s 

where the Hermite operators Fs^^ and the anti-Hermite operators Fs ^ are defined 

by 

Fi±) = i(F,±F]), (3-8) 

-rs=l-15 — , ^ , ^20(+) ' ^21(+) ' ^21(-) ' '°22(+) ' ^22{-) ' 

-^20{+) ' -^2?(+) ' ^2?(-) ' ^22(+) ' ^22(-) ' ^ ^21, -D22}. (3-9) 



Here we use 



^S(±) = ^(^2/4 ± Bi%), {K > 0), (3-10) 



in place of B^j^^ for the quadrupole pairing operators. The interaction strengths Ks 
are given by 

2G?) , 4G(^) , 4G(^) , 4G(^) , 4G(^) , x, 2x, 2x} • (3-11) 

This Hamiltonian is invariant against a rotation by vr about the x-axis. The quantum 
number associated with this is called the signature, r = e~*'^°. The single-particle 
basis with definite signatures are defined by 



4 = -^(cfc + c^), r = -i (a = ^), 

o^fe = 4^(Cfe - Cfc), r = i (a = -J). (3-12) 
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where k denotes the single-particle basis whose magnetic quantum number satisfies 
the condition — 1/2 = even. The operators Fs^^ can be classified according to 
their signature and /C-quantum numbers as shown in Table [H 

The large-amplitude collective vibration responsible for the oblate-prolate shape 
mixing is associated with the K = and 2 components of the interactions in 
the positive-signature (r = +1) sector. Thus the infinitesimal generator of large- 
amplitude collective motion, Q{q), can be written in terms of the single-particle 
basis with definite signature as 

^(«) = E T.'QMil)dldi + QmiMdi, (3-13) 

T kl£T 

where ^ denotes a sum over the signature pairs {k,k), and Q^^^ = Q^. On the 
other hand, the K = 1 component of the interaction in the r = +1 sector and 
the K = 1 and 2 components in the r = — 1 sector contribute to Thouless-Valatin 
equations (|2-25|) . 

Table I. Classification of one-body operators Fi*' in terms of tfie signature r (or a) and K quantum 
numbers. 







r = +l(a = 0) 


r = = 1) 




K = 
K = 
K = 



1 
2 


r4(±) 4(±) o(")(±) r(p)(±) n(+)i 
,-n-p ' -"20(+) ' 20(+) '-^20 J 
rR(n)(±) r(p)(±) n(-)i 
l-°21{-) '-"21{-) '-^21 / 
rR(n)(±) r(p)(±) n( + )i 
t-°22(+) '-"22( + ) '-^22 / 


rR(n){±) r(p)(±) n( + 
l--"21{+) ' -"21(+) ' -^21 
rR(n)(±) r(p){±) n(- 
l-°22(-) '-"22(-) '-'^22 





§4. Solution of the ASCC equations for separable interactions 

4.1. The ASCC equations for separable interactions 

For the separable interactions (j3-7p . the ASCC equations are written aJ^'^^Sli 

H^{q)\hM{q)\<P{q))=^, (4-1) 



5{ct>{q)\ [hM{q\Qm - E4i(9)^i~^ - -B{q)P{q) \m) = 0, (4-2) 

s 



hM{q),-B{q)P{q) 



Y.f'~;:}{q)F^-^^-u\q)Q{q) 



dq 



(4-3) 



where uj'^{q) = B{q)C{q) is the moving-frame QRPA frequency squared, and /ia-/ (g) 
denotes the self-consistent mean-field Hamiltonian in the moving frame, defined by 



hM{q) = Kq)-Y,>^^^HQ)N^^ 



dV_ 
dq 



(4-4) 
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with 

h{q) = ho-Y. {m\ l<^('?)) • (4-5) 

s 

In the above equations, the summation over s is restricted to the operators with 
K = and 2 in the positive-signature sector. We also define the foUowing quantities 

fg'lid) = (<A(g)l [Ft\Q{Q)] \4>{q)) , (4-6a) 
= Ks {Hq)\ [F^^\-B{q)P{q)] \^{q)) , (4-6b) 



.Qiq) 



\Hq)) , (4-6c) 



f(^\q)=B{q)^. (4.6d) 



Note that ah matrix elements are real and (0(g) | Fg ■* \4>{q)) = 0. 

4.2. Overview of solving procedure 

The infinitesimal generators, Q{q) and P{q), which are represented on top of the 
quasiparticle vacuum \(j){q)), are the solutions of the moving-frame QRPA equations, 
while the quasiparticle vacuum \(j){q)), which depends on Q{q), is a solution of the 
moving-frame HB equation. In order to construct the collective path, we have to 
obtain a self-consistent solution for the quasiparticle vacuum and the infinitesimal 
generators. This requires a double iterative procedure for each value of q, because 
the moving-frame HB equation is also solved by iteration. 
Step 0: Starting point 

The shape coexistence phenomena imply that several solutions of the static HB 
equation exist representing different local minima in the potential energy surface. 
We can choose one of the HB solutions and assume that it is on the collective 
path. This starting state is denoted as \4){q = 0)). In the calculation for ^^Se and 
^^Kr in this paper, we choose the HB state at the lowest minimum, which possess 
the oblate shape. As discussed in Ref. '56), a gauge fixing is necessary to solve the 
moving-frame QRPA equations. We choose the "ETOP" gauge. 

Step 1: Initial setting 

Assume that the solution of the ASCC equations at q — Sq is obtained. In order to 
calculate the solution at q, we start from solving the moving-frame HB equation 
(|4Tp . As an infinitesimal generator in the moving-frame Hamiltonian, we use 
an initial trial generator Q{q)^'^^ constructed from the lowest two solutions of the 
moving-frame QRPA equations at q — 6q of the following form 

Q{qf ^ = (1 - e)Qi{q - 5q) + eQ2{q - 5q), (4-7) 

where Qi{q — 6q) and Q2{q — Sq) denote the lowest and the second-lowest solutions 
of the moving-frame QRPA equations at q — 6q, respectively. The mixing param- 
eter e is set to 0.1. This choice is crucial to find a symmetry breaking solution in 
the moving-frame QRPA equations when the moving- frame HB state \(f){q — 5q)) 
and the moving-frame QRPA mode Qi{q — 5q) hold the axial symmetry.^^ The 
5q is set to 0.0157 in the present calculation. 
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Step 2: Solving the moving-frame HB equation 

Using the operator Q^^~^\q) {n > 1), solve the moving-frame HB equation at q 

5(</.W(<7) ^- J^(AM(g))WivW-_ (g)Q(-i)(g) <A(")(g))=0, (4-8) 

T 

with three constraints from the canonical variable conditions 

(0(")(g)|ivW|0(")(g))=O, (4-9) 

if''\q)\Q{q-bq)Y'^\q))=bq. (4-10) 

This step is discussed in ^4.31 in detail. 

Step 3: Solving the moving-frame QRPA equations 

Using the moving-frame HB state \(p^'^\q)), the Lagrange multipliers {X^'^\q))^"'^ 
and dV/dq{q)^"'^ obtained in the previous step, we solve the moving- frame QRPA 
equations with the same gauge-fixing condition as used at the HB state in Step 
0. This determines the infinitesimal generator Q^'^\q) as the lowest solution of 
Eqs. (j4-2p and (j4-3p . Details of this step are described in ^4.41 and Appendix [Bl 

Step 4^ Achieving the self- consistency 

Updating the operator Q^^\q), we go back to Step 2, and repeat Steps 2 and 3 
until all quantities at q converge. 

Step 5: Progression 

Change q io q + 5q and return to Step 1. 
Carrying out Steps 1-5, we obtain a collective path starting from the HB min- 
imum to one direction (q > 0). We then change the sign of 5q and repeat the 
above procedure to the opposite direction (q < 0). In this way we obtain an entire 
collective path. 

After obtaining the solutions of the ASCC equations, we solve the Thouless- 
Valatin eq. (j2-25p at every point on the collective path using the moving-frame HB 
state \4>{q)) to evaluate the rotational moments of inertia Ji{q). Details of this 
calculation are described in Appendix [Cl 

4.3. The moving-frame HB equation in the quasiparticle representation 

The quasiparticle operators a|((g) and 0^(9) associated with the moving- frame 
HB state \(t){q)) are written in terms of the nucleon operators, d\ and d^, with definite 
signature, as 



Its inverse transformation is 



(4-12) 



The U and V matrices are determined by solving the moving-frame HB equation 
(j4-ip . Note that superscripts r(= n,p) for ?7, V , and Fermion operators are omitted 
for simplicity of notation. 
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The moving-frame Hamiltonian is written as 

hM{q) = E E' {ihi\lil)(4dl + idj) - 4?(«)(44 + dldk)) , (4-13) 

T kler 

where the particle-hole part and the particle-particle part of the moving-frame 
Hamiltonian are given by 

{hi\'{q) =h^u){q) - \^^\q)6u' - ^{q)QP{q), (4-14) 

hP =e\^^Su' {4>{q)\ \m {l\F^^¥), (4-15) 

s£ph 

4p= E '^a (0(9)1 ^i+M</>(«))(0|^i+)KrO- (4-16) 

sGpp,hh 

The matrix elements {k\Fs~^^\l) are defined by 

{1\fM\P) = (0|d,F(+)4|0), (0|Fi+)|/rO = (0|F(+)44|0), (4-17) 

where |0) denotes the vacuum for nucleon operators. 
The moving-frame HB equation are thus written 

E' {ihi\i'{q)Ui'k{q) + 4?(?)V;,fc(g)) =4"V,,(g), (4-18) 

Wer 

E' {4,\q)Ui,,{q) + {hi\{q)Vi,k{q)) = - E^Vikiq), (4-19) 
ll'er 

where denotes the quasiparticle energy. These equations are solved under the 
following three constraints: 

(0(g)|iV(«)|</>(g)) = iv("\ (4-20) 

{Hq)\N^M4>{q)) = N^\ (4-21) 
{Hq)\Q{q-Sq)\4>{q)) = Sq. (4-22) 

The Lagrange multipliers, X^"'\q) , X^\q) and dV/dq{q), are determined such that 
these constraints are fulfilled. The expectation values in the moving-frame Hamilto- 
nian are updated using U[k and V/^ thus obtained until self-consistency is achieved. 

In the quasiparticle representation, the moving- frame Hamiltonian hM{q), the 
neutron and proton number operators TV^'^) and the operators Fg"^^ with K = and 
2 in the r = 1 sector are written in the following forms: 

huiq) = {m\ hM{q) \m) + E' iE^,{q)B^^{q) + E-^{q)Bf,-^{q)) , (4-23) 

^^"^ =E'<^(/^^)(4p(9) + ^MP(9)) + E'<^(/^/^)(^^^(«) + 

(4-24) 

= {<Piq)\ m) + Y^F^^]{^^u){Al,{q) + A,,{q)) 

+E'^S(/^'^)(^^-('^) + ^-^^(iy^' (4-25) 

HI/ 
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where 

^Ipiq) = al{q)al{q), A^p(g) = au{q)a^{q), B^y{q) = al{q)a„{q). (4-26) 



Explicit expressions for the matrix elements, N^\n^j^\ and , are given in 
Appendix [Al 

We define the monopole-pairing gaps Z\Q^^(g), the quadrupole-pairing gaps 
A^2K=o 2(^1)^ and the quadrupole deformations D^^^^^iQ) by 



A^o\<l) =G^oHH<i)\ A^^^^^^^ \m) (4-27) 

^2^=0,2(9) = {m\ Di%o,2 l^iQ)) • (4-29) 

4.4. The moving-frame QRPA equations 

The infinitesimal generators, Q{q) and P{q), are represented in the quasiparticle 
representation as 

Qiq) =Q^{q) + Q''{q) 

=Y^Q%mAlM + A,,{q)) + Y^Q^^M{B,M) + B-,,{q)). (4-30) 

P{q) =iY^P,,{q){AlM - A,,[q)). (4-31) 

In the following, we discuss how to get the n-th solution of the moving-frame QRPA 
equations in Step 5 assuming that the (n— l)-th solution Q^'^~^\q) is already known. 
For later convenience, we introduce the following one-body operator 

(4-32) 

with 



P 

"'^fj.p) ■ ■ ^'A,s ~ '^'A,s \^^^^) y^pv 

We can express the matrix elements, Q^p{q) and Pf^p{q), using Eqs. ()4-2p and (j4-3 
as 



- i^i?(^/')(Qfp)^"-'M • (4-33) 
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+ Y.N^^){^,'u')fP{q)\ +54(^^, A')E^A^(^''^'')4i('?)' (4-35) 

r ) s 

where f^j^^iq) = fp^siQ) + fjis^Q)- The metrics ^j, (i = 1 ~ 4) are defined by 

5i(/iP,A')^(-^"''?Wv, 52(Mz^,/iV"0^(A^-i)^,,^v> (4-36) 
g3{nu,fi'u') ={£M'^)^^y-,, g4{fijy,fi'i'') ={£ M'^ £) - 6f,f,>6f^-,, (4-37) 

where M and £ are given by 

.2 



+ (Qj,')^""'^ (e,, - \e, + \E-}j 5,-,^{q), (4-38) 



(4-39) 



The quantities ()4-6p and the canonical variable condition (j2-14p are expressed in 
terms of ((^^(g))^"^ and P{q) as 

= - (0(g)| 10(g)) = -2k.(f(;\ (Q^(<z))(")), (4-40) 



where 



=2^,(Fi+\p(g))-K.(4+],(Q^(g))(")), 

(0(9)1 [iV«,P((z)] \m) = -2z(<\P((z)) = 0, 

{X,Y)^Y.'xifii?)Yifii?). 



(4-41) 
(4-42) 

(4-43) 



Substituting Eqs. (l¥iOD . (I¥ITD . (l¥i2D into Eqs. and (l¥35]l . we derive the 

dispersion equation 



(4-44) 





( 


qQ^PR 






s't' 


qPR,Q 


qPR,PR 


qPR,N 






, qN,Q 


qN,PR 


qN,N 





where the matrix elements of S are given by 



qQ,Q _r)/ rp(—) Tp(-)\ u(~> 



(-) 

A,s ' ^A,s')92 



(4-45a) 
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S?/^=2(Fi;),FW),„ (4.45b) 

S?f=2{Fi-lN^/%, (4.45c) 

+ ^1:3)91 - (4-45d) 

S:^^ =2iFjClFiCX + (<i,it;^),. - ^Jss', (4.45e) 

C^'" =2(i^t^ <').3 + (<i, ^ (4-45f) 

^f;^ =(<,^t').4 - (<'4:i').3, (4-45g) 

5f/^=(<,^i!;').3, (4-45h) 

5f/=(<,<)),3. (4.45i) 

The parentheses in the above matrix elements are defined by 

ix,Y)g^= Y.'x{^x^))gi{^x^)^x'u')Y{^,'u>), (i = i~4). (4.46) 



As we mentioned before, the ASCC equations are invariant under the gauge trans- 
formation associate with number fluctuations. The quantities, 

/Q:i=i,2(9) = -2g(^="'^^ {m\ [i(^="'^)(-), Q(g)] \m) , (4-47) 

and f^j^\q), are transformed by (|2.23p as 

fQ^lUD -/y=i,2(9) - 4a(^="'^)Z^(^='^'^)(g), (4.48) 
ff^'^hl) - a^^=-'^^u;\q). (4.49) 

Thus, we have to fix the gauge in solving the dispersion equation (j4.44p . For both 
neutrons and protons, we choose the "ETOP" gauge^^' 

/fcli(^)=0, /S:i=2('?) = 0- (4-50) 

This gauge fixing condition reduces the dimension of the dispersion equations. We 
can then use the submatrix S' of S, where terms related to the anti-Hermite part 
of the monopole pairing operators, (-^^=12)' dropped. From Eq. (|4.44p . the 
moving- frame QRPA frequency squared u)'^{q) is determined by the condition 

detS'{uJ^{q)) =0. (4.51) 

The lowest uj'^{q) solution (including negative values) are considered as the most 
collective mode at q. Note that we consider imaginary uj{q) solutions as well as real 
ones. Once cij^(g) is determined, f{q), (Q^p(9))*-"^ and P^p{q) are obtained by use 
of the normalization condition 

{m\ mHq))^''\p{q)]\m) = mQ^{'i))^''\p{<i)) = i- (4-52) 
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§5. Requantization of the collective Hamiltonian 

5.1. Requantization and construction of wave functions in the laboratory frame 

Solving the basic equations of the ASCC method and Thouless-Valatin equa- 
tions, we obtain the collective Hamiltonian ()2-24p : we can put the collective mass 
B{q)~^ to unity without losing generality because it just defines the scale of measur- 
ing the length of the collective path. We also put the number fluctuation n to zero. 
Requantization is done simply by replacing the classical variables with the quantum 
operators: 

h^U (5-1) 

I oq 

The Schrodinger equation for the requantized collective Hamiltonian is 

+ Viq) I ^iMkiq, ^) = Ei,k^iMk{q, ^^)- (5-2) 




The collective wave function in the laboratory frame, ^iM,k{q-, is a function of 
the collective coordinate q and three Euler angles 17, and specified by the total 
angular momentum /, its projection M on the laboratory z-axis, and the index k 
distinguishing difl'erent quantum states having the same / and M. Note that the 
three components U of the angular momentum oprerator are defined with respect 
to the principal axes (1, 2, 3) = (x', y', z') associate with the moving-frame HB state 

\m)- 

Using the rotational wave functions ^|^^(i7), we can write the collective wave 
functions in the laboratory frame as 



K=-I 



^lMk{q,^) = '^'lKk{q)\l^^<K{^) (5-3) 



= Yl ^iKkiq) {n\IMK) . (5-4) 
K=0 

Here ^'j^k intrinsic wave functions that represent large-amplitude collective vi- 
brations responsible for the oblate-prolate shape mixing. They are specified, instead 
of M, by the projection K of the angular momentum on the intrinsic z'-axis. We 
assume that the intrinsic states have the positive signature. Then, their K and —K 
components are connected by 

^'iKkiq) = i-Y^'i-Kkiq)- (5-5) 

Accordingly, it is convenient to use new rotational wave functions defined by 

mMK) = -^=^==^l^^{^i,^{f2) + i-Y^y^in)) , (5-6) 
and new vibrational wave functions 
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in place of ^. As the & functions are normalized as 



g^2 

df2&i^Ki^)^M'K'i^) = ^j—y^ii'^mm'Skk', (5-8) 



the normalization of the vibrational wave functions is given by 

nKk{Q)^iKk'{q) = Skk'- (5-9) 



K=0 



5.2. Boundary Conditions 

Multiplying a rotational wave function {Q\IMK) to the Schrodinger equation 
(|5-2p from the left and integrating out the Euler angles J7, we obtain the collective 
Schrodinger equation for large-amplitude vibration: 



^ ^ K'=0 

(5-10) 



where Zot = E^ If /Miq))- 

The boundary conditions can be specified by projecting the collective path to 
the (/?, 7) plane and by using the well-known symmetry properties of the Bohr- 
Mottelson's collective Hamiltonian.'^''^ The deformation parameters P and 7 are 
defined by 

cos 7(g) =x'{m\ Dip \m) /(^o6'), (5-11) 
P{q)smjiq) =V2x' {cP{q)\ D^V \Hq)) /(f^ob^)- (5-12) 

and measure the magnitude and triaxiality of quadrupole deformation of the HB 
mean-field in the moving-frame as functions of the collective coordinate q. Here, 
hujQ denotes the frequency of the harmonic-oscillator potential, x' = X^^i ^-nd the 
harmonic-oscillator length parameter b is related to the radius parameter tq by 

b' = ll^friAl (5.13) 

The boundary conditions for the vibrational collective wave functions depend 
on the character of the collective path in the (/?,7) plane. As we will discuss later 
in ^6.31 the collective path for ^^Se passes through the 7-direction (see Fig(T]). In 
this case, the following boundary conditions are employed at the prolate and oblate 
limits. At the prolate limit, j{qpio) = 0°, the vibrational wave functions are required 
to satisfy 

^iKkiqpvo -q) = {-)^^iKk{qpro + q), (5-14) 

which is equivalent to 

^lKkiqpro) = iK = 2,6,---), (5-15) 
d'^lKk 



dq 



= {K = 0,4:,---). (5-16) 

g='?pro 
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At the oblate limit, 7((7ob) = 60°, the HB mean-field is symmetric about the intrinsic 
y'-axis, and then the boundary conditions are given bj^^ 



(5-17) 

In the case of '^^Kr, the collective path connecting the oblate and prolate shapes 
are not periodic with respect to the 7-direction (see Fig. [3]). Accordingly, we set the 
box boundary conditions at the edge of the path: 

^IKkiqmin) = ^IKkiqmax) = 0. (5-18) 

The matrix elements {IMK\ Trot \IMK') of the rotational kinetic energy oper- 
ator in Eq. (j5-10p can be easily calculated: 

{IMK\ Trot \IMK) =a{q)I{I + 1) + (5-19) 

{IMK\ trot \IM, K + 2)= {IM, K + 2\ f.^t \IMK) 

=c{q){{I + K + 2)(/ + K+ - K)il -K- 

(5-20) 



where 



<;(«) =1 - ^) ■ (5-23) 



.Ml) Ml). 

The other matrix elements are zero. 

5.3. Electric quadrupole moments and transitions 

To evaluate electric quadrupole (E2) moments and transition probabilities, we 
need to derive expressions of the E2 operator in the collective subspace. This can 
be easily achieved by the same procedure as we have derived the quantum collec- 
tive Hamiltonian. As described below, we first take expectation values of the E2 
operators with respect to the moving-frame HB state \4>{q,p)) and then apply the 
canonical quantization procedure. 

In accordance with the quadrupole operators (l3-4j) . we define the E2 operators 
in the model space under consideration as 

^r^=E«e?E<(^04Q, (5-24) 

r kler 

D'!lT=\{b'r + bV)^ (5-25) 

where e^^^ are effective charges. Their expectation values in the collective subspace 
are expanded up to second order in the collective momentum p as 

y(E 



Dlf\q,p) = {cp{q,p)\ ) \4>{q,p)) (5-26) 



D'!.f\^) + \DT^^)P'^ (5-27) 
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where 



D 



l{E2) 



{q)={m\Mf^\<l>{q)) 



D 



//(E2), 



iq)=-{^q)\[[D'if\Qiq)],Qiq)mq))- 



(5-28) 
(5-29) 



The quantities D'l^\q,p) are called collective representations of the E2 operators. 
Note that these are defined in the intrinsic frame associated with the moving-frame 
HB mean-field. We now apply the canonical quantization to them. Then, the 
collective coordinate q and the collective momentum p become quantum operators 
acting on vibrational wave functions ^iKkio.)- We call the requantized E2 operators 
"collective E2 operators" and denote them D^j^ . Thus, the E2 matrix elements 
between two collective vibrational states are evaluated as 



{■l>IKk\ Mf^ \<^IK'k' 



dq^KM { D'l,f\q) - l^D:f^>(q)^ ) <Pjj,,,,{q). 



1 d „//(E2). N d 



2dq 



dq 



(5-30) 



We need to calculate these integrals only for vibrational states which satisfies the 
selection rules of the E2 operators. 

The collective E2 operators are defined in the intrinsic frame, and those 

in the laboratory frame Dj^ ' are obtained by 



(5-31) 



As is well known, i?(E2) values and spectroscopic quadrupole moments Q{Ik) are 
given in terms of reduced matrix elements {Ik\\Dj^ \\Ik) as 



B{E2;Ik^ I'k') = (21 + 1)-^ {Ik\\D^^^^\\l'k') 



(5-32) 



Q{Ik) 



IGtt 



(/, M = I,k\ Z)(^2) I M = I,k) 



IQn f I 2 1 
-I 



; {ik\\D(^'^\\ik). 



(5-33) 



These reduced matrix elements can be evaluated by using the Wigner-Eckart theo- 
rem 

/ 2 



(/, M = I,k\ D^^^^ M = I,k' 

and calculating the left-hand side a;^ 
(/, M = /, k\D^^^'^ M = I, k') 



-I I 



{I,k\\D^''^'>\\l',k'), (5-34) 



V(2/ + 1)(2P + 1) ^ I , 



KK'ti 



V(2/ + l)(2/' + l) ^ {^'jKk\Df^^ WvK' 

KK'fi 

i_k( I 2 /'\ / / 2 I'\ 
I IJ [-K fi K'J ■ 



(-) 



(5-35) 
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In the intrinsic frame, the /i = ±1 components of the collective E2 operator vanish 
and those for the /u = ±2 components are equal. Thus we obtain 



{ik\\b^^'^'^\\i'k') 



+ 



K 2 



-K 



V(2/ + l)(2/' + l)(-/ 

K>0 

I 2 1' 

-K-2 2 K 

12 r 



y'Pi,K+2,k\D'S'^\^r 



Kk' 



(^/xfel D'^^'> \^r,K+2,k' 



/(E2) 



(5-36) 



§6. Results of numerical calculation and discussions 



6.1. Details of numerical calculation 

In numerical calculations, we consider two major shells {N^]^ = 3, 4) for protons 
and neutrons and use the same values for the single-particle energies, the monopole 
pairing strength Gq and the quadrupole particle-hole interaction strength x iii 
Ref. lTSp . The single-particle energies are listed in Table HIl The interaction strengths 
were adjusted to approximately reproduce the pairing gaps and the quadrupole 
deformations obtained by the Skyrme-HFB calculation carried out by Yamagami et 



al.EB These values are gJ,"^ 



G, 



in) 



0.299, Cff ^ = 0.309 and x' 



= 0.320 and x' = xb^ = 0.248 MeV for ^^Se; 
0.255 MeV for ''^Kr. The oscillator frequency 



and the radius parameters are set to hw^ = 41. 2A^/^ MeV and ro = 1.2 fm. For the 
quadrupole pairing strength, we use the self-consistent value derived by Sakamoto 
and Kishimoto,'^ 



G 



(T)self 
2K 



1 / 1 



1 



(a/ 



-1 



(6-1) 



where is the quasiparticle energy evaluated by the BCS approximation at the 
spherical shape. Accordingly, G^2q 



Gself 
21 



Gself 
22 • 



The effective charges e^^^ are written as e^^ 



5epo\ for neutrons and e 



eff 



e + (^Cpoi for protons. For simplicity, we use the same polarization charge Jcpoi = 
0.905e for protons and neutrons, which is chosen to reproduce the experimental 



-B(E2;2J" — > OJ") value^ in Kr. Only this data is available for E2 transitions among 
low-lying states in ^^Se and ^^Kr. This value of (Jcpoi seems a little too large and 
needs further investigation. We take into account the momentum dependent term 
in the collective representation of the E2 operators, Eq. (j5-27p . although numerical 
calculations indicate that it gives only a few percent correction, at most, to the main 
term. 

In the present calculation, we ignore the curvature terms (the fourth, fifth and 
sixth term in Eq. (I4-3P ). to reduce the computational costs. We have checked that 
their contributions are negligible. 

In numerical calculation, careful treatment is necessary for the prolate limit, as 
the moment of inertia about the symmetry axis, J^s^q), vanishes there. Actually, 
this does not cause a problem, because the K components of the vibrational 
wave function also vanish there. To avoid numerical instability, however, we put 
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■J^il) — 10~^^/i^ (MeV)^^, for the prolate limit, and confirmed that this recipe 
works well without losing numerical accuracy. We applied this recipe also for the 
oblate limit where J2{q) vanishes. Actually, the y'-axis component of the vibrational 
wave function also vanishes there, although it is not directly seen in Eq. (j5-10p in 
which the wave functions are decomposed according to the K quantum numbers 
choosing the ^;'-axis as the quantization axis. 

Table 11. Energies in units of MeV of the spherical single-particle levels used in the calculation. 
These values are taken from Ref . '75) . 



orbits 


1/7/2 


2p3/2 1/5/2 


2pi/2 


I59/2 


2d5/2 


lfl7/2 


3Sl/2 


2rf3/2 


neutrons 


-9.02 


-4.93 -2.66 


-2.21 


0.00 


5.27 


6.36 


8.34 


8.80 


protons 


-8.77 


-4.23 -2.41 


-1.50 


0.00 


6.55 


5.90 


10.10 


9.83 



6.2. Properties of local minima in Se and Kr 

We summarize in Table Hill the results of calculation for the properties of the HB 
equilibrium states (local minima in the potential energy surface). For both ^^Se and 
^■^Kr, the lowest HB minimum possesses an oblate shape, while the second minimum 
is prolate. The energy differences between the oblate and prolate minima evaluated 
using the P-l-Q Hamiltonian with (without) the quadrupole pairing interaction are 
300 keV (196 keV) for ^'^Se and 827 keV (626 keV) for ^^Kr. We find no qualitative 
change in the mean-field properties due to the inclusion of the quadrupole pairing 
interaction. 

Table III. The quadrupole deformations and the pairing gaps A[^^ (in MeV) and A^^k (i^i MeV 
fm^) at the HB local minima in ***Se and '^^Kr, calculated with the P-l-Q Hamiltonian including 
the quadrupole pairing interaction. 



(G2 = Gf) 


/3 


7 


^0 


'^0 


^20 


^20 




^22 


"^'Se (oblate) 


0.30 


60° 


1.17 


1.26 


0.08 


0.09 


0.10 


0.11 


®®Se (prolate) 


0.26 


0° 


1.34 


1.40 


0.14 


0.15 








^^Kr (oblate) 


0.35 


60° 


0.92 


1.06 


0.05 


0.06 


0.06 


0.07 


'^^Kr (prolate) 


0.38 


0° 


1.14 


1.27 


0.19 


0.19 









The QRPA collective modes at the oblate and prolate minima can be classified in 
terms of the projections of angular momenta on the symmetry axis, Ky and = K, 
respectively. Table ITVl summarizes the properties the QRPA collective modes at the 
oblate and prolate minima. In ^^Se, the lowest modes are 7-vibrational {Ky or 
Kz=2) and the second lowest modes are /3-vibrational {Ky or Kz=0) both at the 
oblate and the prolate minima. It is seen that the quadrupole pairing interaction 
lowers their excitation energies without changing their ordering. In ''^Kr, the lowest 
QRPA modes at the two minima are both /3-vibrational, if the quadrupole pairing 
interaction is ignored. Note, however, that the Kz = and 2 modes at the prolate 
second minimum are close in energy and their ordering changes when the quadrupole 
pairing interaction is taken into account, whereas the lowest mode at the lowest 
oblate minimum is always /3-vibrational. 

6.3. Collective path connecting the oblate and prolate minima in ^^Se 

We start by solving the basic equations of the ASCC method from the oblate 
minimum {q = 0) and progressively determine the collective path following the al- 
gorithm outlined in §4.2[ Figure [U illustrates the collective path thus obtained by 
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Table IV. The excitation energies uj (in MeV) and the K quantum numbers of the lowest two QRPA 
modes at the oblate and prolate minima in ^*Se and ^'^Kr. The results of the calculation with 
(G2 ~ and without (G2 = 0) including the quadrupole pairing interaction are compared. 

The K quantum numbers here represent Ky or Kz according to the shape (oblate or prolate). 





G2 = 








tj2 = tj2 










LOl 


Kl 


L02 


K2 




Kl 


L02 


K2 


'''*Se (oblate) 


1.555 


2 


2.342 





1.373 


2 


2.131 





^*Se (prolate) 


1.015 


2 


1.915 





0.898 


2 


1.369 





■^^Kr (oblate) 


1.150 





1.909 





1.239 





2.010 


2 


""^Kr (prolate) 


1.606 





1.674 


2 


1.644 


2 


1.714 







0.1 0.2 0.3 0.4 
/3cos7 

Fig. 1. The collective path for ^*Se calculated with the P+Q Hamiltonian including the quadrupole 
pairing interaction. The path is projected onto the (/3, 7) potential energy surface. The dots in 
the figure indicates the HB local minima. The equipotential lines are drawn every 100 keV. 

projecting it onto the (/J, 7) potential energy surface. The path connects the two 
local minima passing through a potential valley lying in the triaxial deformed re- 
gion. The collective path in ^^Se obtained with the P+Q Hamiltonian including 
the quadrupole pairing interaction is very similar to that obtained in Ref. [75|) in 
which its effect was ignored. As solutions of the ASCC equations we obtain various 
quantities; the canonical collective coordinate q, the quadrupole deformations, (3{q) 

and 7((7), the monopole and quadruple pairing gaps A^Q\q) and A^^Ki^l)^ ^^e collec- 
tive potential V{q), the collective mass M{s{q)), the moving-frame QRPA frequency 
squared uj'^{q), and the three rotational moments of inertia Ji{q). These quantities 
are plotted in Fig. [2] as functions of 7(g). The quadrupole deformation /3(g) is al- 
most constant along the collective path while the triaxial deformation 7(g) varies 
and changes the oblate shape to the prolate shape. It is seen that the quadrupole 
pairing interaction slightly increases the f5{q) values for all values of 7(g). 

The collective mass M{s{q)) plotted in Fig. [2] is defined as a function of the 
geometrical length, ds = d/?^ + 13"^ dj'^, in the (/?, 7) plane: 

Misiq)) = M{q)/{{d[5/dqf + f)\d-i/dqf}. (6-2) 

As explained in ^5.11 we can put M{q) = B{qy^ = IMeV"^ here. We have found 
that the quadrupole pairing interaction increases the collective mass. This enhance- 
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Fig. 2. Results of the calculation for ®*Se. The monopole pairing gap Aj (q), the quadrupole 
pairing gaps (q) and A^\q), the collective potential V{q), the collective mass M{s{q)), the 
rotational moments of inertia Jii^q), the lowest two moving-frame QRPA frequencies squared 
01^(5), the axial quadrupole deformation I3{q) and the canonical collective coordinate q are 
plotted as functions of 7(g). Results of the two calculations using the P+Q Hamiltonian with 
(G2 ~ Gl"^") and without (G2 ~ 0) the quadrupole pairing interaction are compared. 

ment takes place almost independent of 7(g), and is mainly due to the decrease of 
dj/dq along the collective path. 

Because the HB mean-field becomes symmetric about the y'- and z'-axes in 
the oblate and prolate limits, respectively, the rotational moment of inertia about 
the y' (2;')-axis vanishes and the other two moments take the same values at the 
oblate (prolate) minimum. Their 7-dependence is similar to that of the irrotational 
moments of inertia. It is found the rotational moments of inertia increase about 
20~30% by the effect of the quadrupole pairing interaction. This enhancement 
as well as that of the inertial functions M{s{q)) is due to the time-odd pair field 
generated by the quadrupole pairing interaction. 

6.4. Collective path connecting the oblate and prolate minima in '^'^Kr 

As for ^^Se, we have determined the collective path for ''^Kr starting from the 
oblate lowest minimum. The collective path projected onto the (/3,7) plane is shown 
in Fig. [3l and various quantities defined along the collective path are plotted in 
Fig. [3] as functions of q. The collective paths calculated with and without including 





. 4. Results of the calculation for '^^Kr. Notations are the same as in Fig. [2l except that the 
quantities are plotted as functions of q along the collective path. The point q = Q corresponds 
to the oblate minimum, while the prolate local minimum is located around q = 3.3 (g = 3.1) for 
the calculation using the P+Q Hamiltonian with (without) the quadrupole pairing interaction. 
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the quadrupole pairing interaction are similar. Because the lowest mode of the 
moving-frame QRPA equations is /3-vibrational around the oblate lowest minimum, 
the path at first goes along the axially symmetric line. Around (/?, 7) = (0.2,60°), 
the character of the lowest mode changes to 7-vibrational, and the path deviates 
from the axially symmetric line. When the collective path reaches the 7 = 0° line, 
the character of the lowest mode again changes to /3-vibrational. Approaching the 
prolate second minimum, the lowest mode changes its character once more to 7- 
vibrational and the collective path deviates from the 7 = 0° line. 

We note that the lowest two modes at the prolate second minimum are very close 
and their ordering in energy may be sensitive to the interactions used. We examined 
that, for example, the lowest mode at the prolate second minimum becomes (5- 
vibrational, if the quadrupole pairing interaction is switched off, and in this case the 
axial symmetry breaking takes place at a larger /5 value beyond the prolate second 
minimum. In such a situation, two collective coordinates may be needed to describe 
the collective dynamics in a better way. This serves as an interesting subject for 
future investigation. It should be emphasized that such a problem arises only locally 
in a small region on the 7) plane and the collective path is well defined globally. 

At large f3 region beyond the oblate minima {q < 0) along the 7 = 60° line, the 
lowest K = mode exhibits a strong mixture of /3-vibration (fluctuation of axially 
symetric shape) and neutron pairing vibration (fluctuation of pairing gaps), and the 
calculation to find the collective path eventually stops when the neutron monopole 
pairing collapses. 

As for ^^Se, we have found that the collective mass and the rotational moments 
of inertia increase also for ''^Kr due to the time-odd pair field generated by the 
quadrupole pairing interaction. We note that the collective mass M(s{q)) diverges 
at large deformation. This behavior, seen also in the previous work,'^''^''^ is 
associated with the disappearance of the pairing gaps. 

6.5. Excitation spectra and quadrupole transitions in ^^Se 

The collective Schrodinger equation (|5-2|) is solved with the boundary conditions 
(j5-14p and (j5-17p for ^^Se to get energy spectra, quadrupole moments and transi- 
tion probabilities. The result of calculation is displayed in Fig. [5l The calculation 
yields the excited prolate rotational bands as well as the oblate ground state band. 
It is seen that the inter-band E2 transitions are weaker than the intra-band E2 
transitions, indicating that the oblate-prolate shape coexistence picture holds. The 
calculation suggests the existence of the excited 0"*" state which has not yet been 
found in experiment. The spectroscopic quadrupole moments presented in Fig. [6] 
are also consistent with the oblate-prolate shape coexistence picture: The yrast 
states possess positive spectroscopic quadrupole moments indicating the oblate de- 
formation, while the second lowest states for each angular momentum have negative 
value indicating the prolate deformation. In Fig. 0, the excitation spectra calculated 
with and without including the quadrupole pairing interaction are compared. We 
see that the quadrupole pairing plays an important role in decreasing the excitation 
energies. This is because the time-odd pair field generated by the quadrupole pairing 
enhances the collective mass and the rotational moments of inertia. 

In Fig. [71 the vibrational wave functions are presented. One may notice that 
the behaviors of the O"*" states are significantly different from the 1^0 states: The 
vibrational wave functions of the lowest and the second lowest 0"*" states spread over 
the entire collective path indicating that the oblate and prolate shapes are strongly 
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Fig. 5. Excitation spectra and _B(E2) values of low-lying states in ®*Se calculated by means of the 
ASCC method. In the left (middle) panel, the quadrupole pairing is ignored (included) in the 
microscopic Hamiltonian. Experimental datcP are displayed in the right panel. The _B(E2) 
values larger than 1 W.u. are indicated in parentheses besides the arrows in units of fm*. 
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Fig. 6. Spectroscopic quadrupole moments of low-lying states in ^**Se. The left and right panels 
show the spectroscopic quadrupole moments of the yrast states and of the second lowest states 
in each angular momentum, respectively. The units for the right panels are indicated besides 
the right vertical lines. Results of calculation with (without) including the quadrupole pairing 
interaction are indicated by filled (open) squares. 
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Fig. 7. Vibrational wave functions 'l>iKk{<l) of the yrast states {left) and the second lowest states 
in each angular momentum {right): in ^*Se, In each panel, different K components of the 
vibrational wave functions are plotted as functions of 7(g). The calculation is performed with 
the P+Q Hamiltonian including the quadrupole pairing interaction. 



mixed via the triaxial degree of freedom. In contrast to the 0"'' states, the 1^0 
wave functions contain the K ^ components which take the maximum values at 
the oblate limit. We can see this trend more clearly by plotting the collective wave 
functions squared. This is presented in Fig. [8j The vibrational wave function of the 
ground 0"*" state spreads over the entire region of 7, while that of the excited 0"*" 
state exhibits prominent peaks both at the oblate and prolate limits. In contrast, 
the vibrational wave functions of the 1^0 yrast states are localized around the 
oblate shape, while those of the second lowest states (in each angular momentum) 
are localized around the prolate shape. The localization develops with increasing 
angular momentum. In the yrast states, all the K ^ components exhibit the 
maxima at the oblate shape, while the K = component dominates at the prolate 
shape in the second lowest states. 

In order to evaluate the oblate-prolate shape mixing in a more quantitative 
manner, we define the oblate and prolate probabilities as follows: 

rqo ^ /-qmax ^ 

lKkiQ)\'^, (6-3) 

where we assume that gmin < ^ob < lo < Qpm < ^max- The "boundary" between the 
oblate and the prolate regions is set to the top of the potential barrier between the 
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Fig. 8. Vibrational wave functions squared |^/Kfe(g)P the yrast states {lef£) and the second lowest 
states in each angular momentum {right): in ^*Se, In each panel, different K components are 
plotted as functions of 7(g). The solid (dashed) line shows the result of calculation using the 
P+Q Hamiltonian with (without) the quadrupole pairing interaction. Note that the vibrational 
wave functions are normalized as Eq. (|5-9|) with respect to the collective coordinate q, so that 
the factor d'y/dq is multiplied to them when integrating with respect the triaxial deformation 
parameter 7. The d'y/dq values calculated including the quadrupole pairing interaction are 
larger than those without including it for all range of q. 

two minima, or at 7 = 30°. Figure [9] shows these probabihties for ^^Se. The oblate 
and prolate states are strongly mixed in the O"*" states. It is clearly seen that the 
shape mixing rapidly decreases as the angular momentum increases. 

6.6. Excitation spectra and quadrupole transitions in '^^Kr 

For ^^Kr, the collective Schrodinger equation is solved under the boundary con- 
ditions (|5-18p . The result of calculation exhibits two coexisting rotational bands; 
see the energy spectra and the B(E2) values displayed in Fig.[10l The spectroscopic 
quadrupole moments presented in Fig. [11] indicate that the yrast band possesses the 
oblate character, while the excited band has the prolate character. For all states 
including the 0"*" states, the inter-band B(E2) values are smaller about one-order 
of magnitude than the intra-band B(E2) values, and they rapidly decrease as the 
angular momentum increases. This indicates that the oblate-prolate shape mixing 
is rather weak in '^^Kr. 

In Fig. [T2\ the vibrational wave functions are plotted. It is seen that the wave 
function of the 0^ state is well localized in the oblate region, while that of the O2 



28 



N. Hinohara, T. Nakatsukasa, M. Matsuo and K. Matsuyanagi 





1 


Gi=Q 














1 






0.8 












• 


• 


0.8 






U.D 
















n R 

U.D 




P(l,k) 


0.2 






■ 


f 


--■ 






'""G 


o 


0.2 



P{l,k) 




1 
















1 






0.8 
0.6 








1 


■ 


• 


• 

--■ 


0.8 
0.6 




P(l.k) 


0.4 


• 












0.4 


P(l,k) 




0.2 


















0.2 









2 


4 


6 





2 


4 


6 







/(g.s.) /(1st) 

Fig. 9. The oblate and prolate probabilities evaluated for individual eigen-states in ®*Se. The 
upper (lower) panel shows the probabilities calculated using the P+Q Hamiltonian without 
(with) the quadrupole pairing interaction. The probabilities defined by setting the boundary 
at the barrier top (7 — 30°) are shown by squares (circles). 

state exhibits the major peak in the prolate region. In the yrast states with / 7^ 0, 
the locahzation about the oblate shape further develops for all X-components of 
the vibrational wave functions. The magnitude is larger for higher K. In contrast, 
the collective wave functions of the second lowest states in each angular momentum 
are essentially composed of the X = component which localize in the prolate 
region. Figure [13] shows the vibrational wave function squared. Rather weak oblate- 
prolate shape mixing is seen only for the excited 0"*" state, and other members of the 
rotational bands possess well-defined oblate or prolate characters. The oblate and 
prolate probabilities are presented in Fig. [TH The shape mixing in the 0"*" states is 
much weaker compared to ^^Se, and almost diminishes at finite angular momentum. 

§7. Concluding Remarks 

The shape coexistence/mixing phenomena in the low-lying states of ^^Se and 
^^Kr have been investigated using the ASCC method. The excitation spectra, the 
spectroscopic quadrupole moments and the E2 transition properties of the low-lying 
states have been evaluated for the first time using the ASCC method. We have 
derived the quantum collective Hamiltonian that describes the coupled collective 
motion of the large-amplitude vibration responsible for the oblate-prolate shape 
mixing and the three-dimensional rotation of the triaxial shape. The calculation has 
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Fig. 10. Same as Fig. [5] but for ^^Kr. Experimental data are taken from Refs. [64 |l .165 11 . iGGjl 
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Fig. 11. Same as Fig. [6] but for ^^Kr. See the caption of Fig. [6] 



yielded the excited prolate rotational band as well as the oblate ground-state band. 
The basic pattern of the shape coexistence/mixing phenomena has been well re- 
produced using the one-dimensional collective path running on the two-dimensional 
(/9, 7) plane. This collective path is self-consistently extracted from the huge dimen- 
sional TDHB manifold. Thus, the result of calculation indicates that the TDHB 
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Fig. 12. Vibrational wave functions <PiKk{q) of low- lying states in ^^Kr plotted as functions of q. 
See the caption to Fig. [T] 



collective dynamics of the shape coexistence/mixing phenomena in these nuclei is 
essentially controlled by the single collective coordinate microscopically derived by 
means of the ASCC method. 

We have also shown that the low-lying states are significantly better described 
by including the quadrupole pairing interaction. The reason is that the time-odd 
component of the mean-field generated by the quadrupole pairing interaction en- 
hances the collective mass of the vibrational motion and the moments of inertia of 
the rotational motion, lowering the energy of the collective excitation. 

The present calculation clearly indicates that the oblate-prolate shape mixing 
decreases as the angular momentum increases. This implies that the rotational 
dynamics plays the major role to realize the localization of vibrational wave functions 
around the oblate and prolate minima under the situation that the barrier between 
these local minima is very low. We shall attempt a more detailed investigation on 
the dynamical reason why the rotational motion hinders the oblate-prolate shape 
mixing in a separate paper. 
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Fig. 13. Vibrational wave functions squared \$iKk{q)\^ of low-lying states in '^^Kr plotted as func- 
tions of q. See the caption to Fig. [8l 



Appendix A 

Quasiparticle Representation of One-Body Operators 



Since the moving mean- field \(j){q)) has the positive signature, the conditions 

hold. The matrix elements of the pairing one-body operators F^^^ 2 3 6 8 ii '^ith 
K = {),2 and r = +1 {A^^^^^\b[^^I^^^ , and ^^2(+?) ™ ^q. (I¥25]) are 



kl 
kl 

=Y^{^\F^^^k^){VJ^{q)UkM) ± UkMyM)- 



(A-2) 
(A-3) 
(A-4) 

The matrix elements of the particle- hole operators with = 0, 2 and r = {D^o 
and b''^^) in Eq. (I¥25]l are 



kl 



(A-5) 



kl 



32 



N. Hinohara, T. Nakatsukasa, M. Matsuo and K. Matsuyanagi 



P(l,k) 



1 

0.8 « 

G2=0 

0.6 
0.4 

0.2 , 
■ 

1 

0.8 



0- 



/-tself 
L/2=tJ2 



P(l.k) 



0.6 
0.4 
0.2 





2 4 
/(g.s.) 



6 



■ 





1 

0.8 

0.6 

0.4 

0.2 



1 

0.8 
0.6 
0.4 
0.2 




P(l.k) 



P(l.k) 



2 4 6 
/(1st) 



Fig. 14. The same as Fig. [9] but for ^^Kr. 



A,s=13,15 



kl 



kl 

The matrix elements of the particle number operators in Eq. (j4-24p are 



N(-\q) = {<P{q)\ iVW 10(g)) = 2^'^%{qf 

k^T A 

^(^^) =Y^{Uk^{<i)ykv{q) - uj,My-kM)^ 

k&T 

N^i^^^) =Y^{UkMUku{q) - Vj,M)y-kM)- 



(A-8) 
(A-9) 
(A-IO) 



k& 



The constraint about the Q{q — Sq) operator (I4-22P in the moving-frame HB equation 
is written as 

{4>{q)\ Q{q - 5q) |0(g)) = Y^Qkl{q - 5q)Y^Vk-MVij,{q). (A-11) 

kl ^l 
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Appendix B 

Determination of the B-part of Q{q) 

We show that the i?-part of the operator {Q{q)) can be determined through 
its A-paxt {Q^{q)), which is obtained by solving the moving- frame QRPA equa- 
tions. In terms of the quasiparticle operators, a| and Ui, defined by the Bogohubov 
transformation 



(B-1) 



the Hermitian operator Q{q) is written as 

= E {qUM'^] + Q^*iQ)'^J^^ + , (B-2) 



where 



= U^QV, = U^QU - V^QV. (B-3) 

Thus the matrices, Q and , can be written in terms of as 

Q =(C7t)-iQ^y-i, 

=Q^V-^U -V^U^y^Q^. (B-5) 

We cannot directly use these relations, however, for determining Q and , because 
the matrices Q and calculated by ()B-4p and ()B-5P are not Hermitian. We have 
to construct Hermitian matrices Q and from Q^. This is achieved by adding a 
symmetric matrix S to the solution of the moving-frame QRPA equation, which we 
here denote Qq, as 

Q^ = Q^ + S. (B-6) 

For this matrix, the operator Q is written as 

Q =(C/t)-i(Q^ + S)V-^ = Qo + {U^)-'SV-\ (B-7) 
Qt =(yt)-i(QA ^ ^)t^-i ^ gt ^ (yt)-i5t[/-i. (B.8) 

From the Hermicity condition, Q = , we obtain the following equation for 5. 

{V^r'S^U-' - [U^r'SV-^ = Qo- qI (B-9) 
which determine the symmetric matrix S. Let us write the above equation explicitly, 
Y,{{V-\,{U~% - {U-\,{V-%}Ski ={Qoh - (.Qo)ji, (B-io) 

kl 

where we assume that all quantities are real. Both the number of unknown quantity 
and number of equations are the same, A^(A^ -|- l)/2, N being the dimension of the 
matrix. Therefore it is possible to determine the matrix S by solving this equation. 
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In the case of the P+Q model, we start from the skew symmetric matrix Q^, 

iQo)MQ) = IqU^i)^ iQoh.iQ) = -^QtM- (B-ii) 

The symmetric matrix S^jp = Sp^ is determined by solving the following equation: 

- (.U~^)^.k{V-%l}S^J, ={Qo)ki - {Qo)ik (B-12) 

E'{(^~')Mfc(f^'')pr- =(Qo)fcr- (Qok (B-13) 



where 



{QM<i) = Y^{u-^)AQt)A'iW%h (B-14) 

m-kM = Y^{U-%j,{Q^),^{q){V-^)^j. (B-15) 



As the relation {Qq)i_ip = {Qq)up. holds, the matrix S satisfies the relation S^p = 
-S^fi, and then, the Eqs. ([B42ll and (lB43]l are written as 

J2'i-(^'^)-k{U-^)^i - {U-\k{V-^)fii}S^,p ={Qo)m - {Qo)ik. (B-16) 

Using the transformed matrix Q^' (q) 

Qi'M = iQo)MQ) + S^.u, (B-17) 
the Hermite matrices Q{q) and Q^{q) are obtained as follows: 

Qkliq) = Y^{U~^)^kQ%{q){V-%u (B-18) 

Q-kM = Y^{u-')-,iQtmy'')uh (B-19) 

QUi) = Y^Uki^Qki{q)Uiu - V-k^Q-^i{q)V,i, (B.20) 

kl 

QUi) = Yl'Up-kQ-kTiQ)UT, - Vk-^QMVip. (B-21) 

kl 

Appendix C 

Calculation of the Rotational Moments of Inertia 

For the separable interactions (j3-7p . the Thouless-Valatin equations (j2-26p de- 
termining the three rotational moments of inertia Ji{q) about the principal axes at 
a non-equilibrium state \4>{q)) can be written in the following form: 



s s 

1 



^J-\q)I,\cPiq))=0, (C-1) 
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where 

&9) = i^s (</>(g)| [F^'-^Mq)] m) , (C-2) 

41(9) = {m\ [Fs^-\nQ)] m) ■ (c-3) 

The quasiparticle representation of the angular momentum operators are 

4 lA,xil^'')i'^lp + ^M^^) + ^'lB,x{HiB^,u - Bfip), (C-4) 

jJLV 

+ Y^lB,y{lJiu){B^^-B-^,), (C-5) 

/if 

+ Y^lB,zil^^)iB^v + Bfi,), (C-6) 

HI/ 

where the matrix elements of Ia,x, lA,y and Ia,z are given by 

kl 

lA,yif^i^) =lA,yifi^) = Y^ik\Iy\l)Uk^iq)VM, (C-8) 

kl 

Ia,z{i^^) = - Ia,z{i^v) = mkUk^{q)V^^{q). (C-9) 

fe 

The residual interactions with (r = +1, iiT = 1), (r = — 1, = 1) and (r = —1,K = 
2) contribute to rotations about the x, y and z-axis, respectively. The quasiparticle 
representation of one-body operators having these quantum numbers are 

=^F^^^{f^u){Al±A,,) 

+ Y^f'^b}{H{B^.-Bj,,), {r = +l,K = l) (C-10) 

+ (r = -l,K = l) (C-11) 

=^Vif,V^){4^ ± - (4p ± ^Ap)} 

+ + {r = -l,K = 2). (C-12) 
The matrix elements of the quadrupole pairing operators are 

^£U9(/^^) ='2Y^mSl-^\mviM^kM ± Uk.mj-Mi (c-i3) 



36 N. Hinohara, T. Nakatsukasa, M. Matsuo and K. Matsuyanagi 

kler 

^£Ui2(/"^) =E'(0|^22(-?l^0(-^^(^)^fc.(9) ± Uk^{q)UUq)), (C-15) 

kler 

Note that the quadrupole operators D22 ^ does not contribute to the moments 

of inertia. 

The three Thouless-Valatin equations, (jC-ip . at non-equihbrium can be solved 
independently. The angle operators, ^xiq),'I'y{q) and 'I^z{q), can be written as 

Mq) = iY^'I'A,x{f^^)iAlp - ^t^p) + {B- part), (C-17a) 

^yil) = Y^'^A,y{^J^^^){Al^ + A^^) + ^A,y{l^v){Al^ + Aj,^) + (S-part), (C-17b) 

Mq) = ^^'-^A.^lHl^i^ - A^^) + a^^,^(Ai^)(Atp - App) + (S-part). (C-17c) 
These matrix elements are easily obtained from Eq. (IC-ip as 

(- E - J.~Hq)lAAH^ , 

[- E - J.-\q)iAA-^iy)^ ■ 

It is easy to confirm that f!^Jg{q) = 0, = and f^Jg{q) = 0. Substitut- 
ing (jC-lSP into the quantities, f!^\{q),f!^ ^s(l) ^ f^^^il) ^ ^^'^ canonical variable 
condition (j2-26p . we obtain 

= 2ks E(^iti' ^A!i)i^+4S('?) + '2'^s(.F[^llA,.)E+J-Hq), (C-19a) 

s' 

ftJsi<i) = -^s{m\[F^-\^ymq)) 

=2Ks^{F^-lF[-],)E-i-jAQ)-^'^siF['^ (C-19b) 

it^q) =i^s {Hq)\[F^^\Mq)m.q)) 

=2ks ^{F[^],TZ)E-i^l,{q) + 2Ks(,F[+],lZ)E-J,^Hq), (C-19c) 



'FA,y{fJ-A 
'FA,y{fii^) 

'FAAf^'') 
'FaAI^'^) 





-1 


F, 


+ F^ 




1 


F, 


+ F^ 




1 


F, 


+ F^ 




1 


F, 


+ F^ 




1 


F-, 





(C-18a) 
(C-18b) 
(C-18c) 
(C-18d) 
(C-18e) 
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s 



= 1, (C-20a) 



= 2Y,{F\J,lA,y)E-r^jA<l) - 2{lA,y,lA,y)E-JyHq) 
s 

= -1, (C-20b) 

fiu flu 

= 2E(^i!i'^)^-/S'(9) + 2{lA,z,7Z)E-J;'Hq) 

s 

= 1, (C-20c) 



where 



= X{nu) - X{up), (C-22) 



Equations ()C-19P and ()C-20p are Unear equations with respect to f!^l{q) and ^{q), 
and can be rewritten as follows: 

^ / 2^,(fW, - 6,, 2kAf[':IIa,x)e+\ ( i^iM)] = f 0\ 

^\ 2{F^+],Ia,.)e+ 2{Ia,.,Ia,x)e+, J [jc^Hq) J Vv' 

(C-24a) 

2^.(^1^ , - Sss' -2KAFi-]jZy)E^\ ( a-jm] = (0\ 

-2{F[-]jA,y)E- 2{lA,y,lA,y)E~ )\^y' I W' 

(C-24b) 



E 



2(Fi+U^,,)£;_ 2{Ia,.,Ia,.)e- J \Jz-\q)J \^>' 



(C-24c) 
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